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Abstract 

In situ measurements of ice crystal concentrations and sizes made with aircraft instrumentation 
over the past two decades have often indicated the presence of numerous relatively small (< 50 pm 
diameter) crystals in cirrus clouds. Further, these measurements frequently indicate that small 
crystals account for a large fraction of the extinction in cirrus clouds. The fact that the instruments 
used to make these measurements, such as the Forward Scattering Spectrometer Probe (FSSP) and 
the Cloud Aerosol Spectrometer (CAS), ingest ice crystals into the sample volume through inlets has 
led to suspicion that the indications of numerous small-crystals could be artifacts of large-crystal 
shattering on the instrument inlets. We present new aircraft measurements in anvil cirrus sampled 
during the Tropical Composition, Cloud, and Climate Coupling (TC4) campaign with the 2- 
Dimensional Stereo (2D-S) probe, which detects particles as small as 10 pm. The 2D-S has detector 
"arms" instead of an inlet tube. Since the 2D-S probe surfaces are much further from the sample 
volume than is the case for the instruments with inlets, it is expected that 2D-S will be less 
susceptible to shattering artifacts. In addition, particle inter-arrival times are used to identify and 
remove shattering artifacts that occur even with the 2D-S probe. The number of shattering artifacts 
identified by the 2D-S interarrival time analysis ranges from a negligible contribution to an order of 
magnitude or more enhancement in apparent ice concentration over the natural ice concentration, 
depending on the abundance of large crystals and the natural small-crystal concentration. The 2D-S 
measurements in tropical anvil cirrus suggest that natural small-crystal concentrations are typically 
one to two orders of magnitude lower than those inferred from CAS. The strong correlation between 
the CAS/2D-S ratio of small-crystal concentrations and large-crystal concentration suggests that the 
discrepancy is likely caused by shattering of large crystals on the CAS inlet. We argue that past 
measurements with CAS in cirrus with large crystals present may contain errors due to crystal 
shattering, and past conclusions derived from these measurements may need to be revisited. Further, 
we present correlations between CAS spurious concentration and 2D-S large-crystal mass from 
spatially uniform anvil cirrus sampling periods as an approximate guide for estimating quantitative 
impact of large-crystal shattering on CAS concentrations in previous datasets. We use radiative 
transfer calculations to demonstrate that in the maritime anvil cirrus sampled during TC4, small 
crystals indicated by 2D-S contribute relatively little cloud extinction, radiative forcing, or radiative 
heating in the anvils, regardless of anvil age or vertical location in the clouds. While 2D-S ice 
concentrations in fresh anvil cirrus may often exceed 1 cm -3 , and are observed to exceed 10 cm -3 in 
turrets, they are typically ~0.1 cm -3 and rarely exceed 1 cm -3 (<1.4% of the time) in aged anvil cirrus. 
We hypothesize that isolated occurrences of higher ice concentrations in aged anvil cirrus may be 
caused by ice nucleation driven by either small-scale convection or gravity waves. It appears that the 
numerous small crystals detrained from convective updrafts do not persist in the anvil cirrus 
sampled during TC-4. 
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Abstract. In situ measurements of ice crystal concentra- 
tions and sizes made with aircraft instrumentation over the 
past two decades have often indicated the presence of nu- 
merous relatively small ( < 50 /xm diameter) crystals in cirrus 
clouds. Further, these measurements frequently indicate that 
small crystals account for a large fraction of the extinction 
in cirrus clouds. The fact that the instruments used to make 
these measurements, such as the Forward Scattering Spec- 
trometer Probe (FSSP) and the Cloud Aerosol Spectrometer 
(CAS), ingest ice crystals into the sample volume through 
inlets has led to suspicion that the indications of numerous 
small-crystals could be artifacts of large-crystal shattering 
on the instrument inlets. We present new aircraft measure- 
ments in anvil cirrus sampled during the Tropical Composi- 
tion, Cloud, and Climate Coupling (TC4) campaign with the 
2-Dimensional Stereo (2D-S) probe, which detects particles 
as small as 10 /xm. The 2D-S has detector “arms” instead 
of an inlet tube. Since the 2D-S probe surfaces are much 
further from the sample volume than is the case for the in- 
struments with inlets, it is expected that 2D-S will be less 
susceptible to shattering artifacts. In addition, particle inter- 
arrival times are used to identify and remove shattering ar- 
tifacts that occur even with the 2D-S probe. The number of 
shattering artifacts identified by the 2D-S interarrival time 
analysis ranges from a negligible contribution to an order 
of magnitude or more enhancement in apparent ice concen- 
tration over the natural ice concentration, depending on the 
abundance of large crystals and the natural small-crystal con- 
centration. The 2D-S measurements in tropical anvil cirrus 
suggest that natural small-crystal concentrations are typically 
one to two orders of magnitude lower than those inferred 
from CAS. The strong correlation between the CAS/2D-S ra- 
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tio of small-crystal concentrations and large-crystal concen- 
tration suggests that the discrepancy is likely caused by shat- 
tering of large crystals on the CAS inlet. We argue that past 
measurements with CAS in cirrus with large crystals present 
may contain errors due to crystal shattering, and past conclu- 
sions derived from these measurements may need to be re- 
visited. Further, we present correlations between CAS spu- 
rious concentration and 2D-S large-crystal mass from spa- 
tially uniform anvil cirrus sampling periods as an approxi- 
mate guide for estimating quantitative impact of large-crystal 
shattering on CAS concentrations in previous datasets. We 
use radiative transfer calculations to demonstrate that in the 
maritime anvil cirrus sampled during TC4, small crystals in- 
dicated by 2D-S contribute relatively little cloud extinction, 
radiative forcing, or radiative heating in the anvils, regardless 
of anvil age or vertical location in the clouds. While 2D-S ice 
concentrations in fresh anvil cirrus may often exceed 1 cm -3 , 
and are observed to exceed 10 cm -3 in turrets, they are typ- 
ically ~0.1 cm -3 and rarely exceed 1 cm -3 (<1.4% of the 
time) in aged anvil cirrus. We hypothesize that isolated oc- 
currences of higher ice concentrations in aged anvil cirrus 
may be caused by ice nucleation driven by either small-scale 
convection or gravity waves. It appears that the numerous 
small crystals detrained from convective updrafts do not per- 
sist in the anvil cirrus sampled during TC-4. 


1 Introduction 

Accurate information about cirrus ice crystal size distribu- 

tions is important for a number of reasons. The balance be- 
tween cirrus longwave and solar cloud radiative forcing de- 
pends, among other factors, on ice crystal size (Stackhouse 

and Stephens, 1991; Jensen et al., 1994). Cloud remote- 

sensing measurements typically measure only part of the size 
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distribution. In order to retrieve parameters such as ice wa- 
ter content, apriori information about the shape of ice crys- 
tal size distributions is required. Sanderson et al. (2008) re- 
cently used an ensemble of thousands of general circulation 
model (GCM) simulations with perturbed cloud parameteri- 
zations to evaluate the relative impact of different cloud pa- 
rameters on climate sensitivity in the model. One of the two 
leading factors was ice crystal fallspeed, which is a function 
of the ice crystal size assumed in the cirrus parameteriza- 
tions. Mitchell et al. (2008) showed that cirrus ice mass, cir- 
rus coverage, and cloud radiative forcing simulated in a cli- 
mate model are very sensitive to assumptions about the con- 
centration of small (diameter < 60 /x m) crystals in the cirrus 
parameterization. Lastly, if we are to evaluate potential im- 
pacts of anthropogenic aerosols on cirrus microphysical and 
radiative properties, we must first have reliable information 
about ice crystal sizes and concentrations. 

Airborne in situ instruments have been used on aircraft to 
measure the sizes and concentrations of ice crystals in cir- 
rus clouds for more than three decades. If we focus on the 
small end of the cirrus ice size spectrum ( <50 /xm diame- 
ter), the vast majority of the measurements were made with 
FSSP-type (Forward-Scattering Spectrometer Probe) instru- 
ments that measure the scattering of light from the crys- 
tals, and retrieve the crystal size using Mie scattering the- 
ory. Instruments based on this technique include the FSSP- 
100, FSSP-300, SPP-100, MASP (Multi- Angle Spectrome- 
ter Probe), and the CAS (Cloud Aerosol Spectrometer) com- 
ponent of CAPS (Cloud, Aerosol, Precipitation Spectrome- 
ter) (Baumgardner et al., 2001). More recently. Counterflow 
Virtual Impactors (CVIs) have been used to infer ice crystal 
concentrations (Strom et al., 1997) by counting the residual 
particles after ice sublimation and assuming one residual par- 
ticle per ice crystal. These in situ measurements have been 
used to characterize the microphysical and radiative proper- 
ties of both midlatitude cirrus (e.g. Ivanova et al., 2001; Law- 
son et al., 2006a) and tropical cirrus (e.g. Platt et al., 1989; 
Knollenberg et al., 1993; Garrett et al., 2005). The com- 
mon aspect of the above instruments is that they all ingest 
ice crystals through an inlet and are flown on jet aircraft (for 
cirrus measurements) at high speeds, raising the possibility 
of small-particle artifacts due to large-crystal shattering. 

The measurements of abundant small crystals have been 
questioned for several reasons. As an example, the CAS 
measurements in anvil cirrus from the Cirrus Regional Study 
of Tropical Anvils and Cirrus Layers - Florida Area Cirrus 
Experiment (CRYSTAL-FACE) are suspicious for a num- 
ber of reasons: first, abundant small crystals are indicated 
throughout the clouds (in fresh anvils, in aged anvils, and 
throughout the vertical extent of the clouds) (Garrett et al., 
2005). Size sorting as anvils age should produce regions in 
the lower parts of the clouds where large crystals are present 
with few small crystals. Second, effective radii derived from 
the CAPS measurements are typically in the 5-20 /im range; 
these values are substantially smaller than those derived from 


remote-sensing measurements (typically 30-40 /im) (Hong 
et al., 2007). Lastly, there is a strong correlation between 
small-crystal concentration derived from CAS and large- 
crystal mass derived from imaging probes such as CIP or 
2D-C (Heymsfield, 2007). Such a correlation to be expected 
if shattering of large crystals dominates the measured con- 
centration of small crystals. 

We refer here to crystals with lengths less than 50 /x m as 
“small crystals”. Past studies debated the role of small crys- 
tals for tropical anvil cirrus albedo, defining “small” as less 
than 20-25 /x m (Knollenberg et al., 1993; Zender and Kiehl, 
1994; Heymsfield and McFarquhar, 1996; McFarquhar and 
Heymsfield, 1997). Some of these studies indicated even 
these very small crystals were important, while others con- 
cluded the small crystals had little impact on cirrus albedo. 
We use a larger upper limit (50 /xm) here because the entire 
size range from about 10 /xm to the upper size limit of the 
CAS spectrometer seems to be affected by crystal shattering 
(McFarquhar et al., 2007). 

Past studies have provided mixed evidence about the plau- 
sibility of large concentrations of small crystals (>lcm -3 
crystals smaller than ~50 /xm) in ice clouds. Gardiner and 
Hallett (1985) compared FSSP measurements with replica- 
tor measurements and concluded that the FSSP-measured ice 
concentrations could be erroneously enhanced by orders of 
magnitude in the presence of large, irregular ice crystals. 
However, Ivanova et al. (2001) argued that the counting ef- 
ficiency of the replicator decreases with decreasing crystal 
size such that one would expect undercounting of small crys- 
tals. Gayet et al. (1996) compared FSSP and 2D-C mea- 
surements and concluded that FSSP concentrations are ex- 
aggerated in cirrus with large crystals. However, Gayet et al. 
(2002) and Lawson et al. (2006a) measured FSSP concen- 
trations of small ice greater than 1 cm -3 in midlatitude cir- 
rus in the absence of crystals larger than 100 /xm. Field 
et al. (2003) evaluated the interarrival times for crystals in the 
FSSP instrument sampling midlatitude cirrus and concluded 
that the crystal concentrations can be exaggerated by factors 
of 2-3 and occasionally by as much as a factor of 5. We 
note that much larger crystals were present in the anvil cir- 
rus discussed here than in the midlatitude cirrus analyzed by 
Field et al. (2003). As shown by Korolev and Isaac (2005), 
the number of fragments produced by shattering increases 
rapidly with the size of the large crystals. Hence, the ex- 
aggeration of small-crystal concentration due to shattering 
is likely to be much larger for the tropical anvil cirrus dis- 
cussed here than the Field et al. (2003) analysis suggested 
for midlatitude cirrus. Heymsfield (2007) recently compared 
extinction and ice water content measured with FSSP and 
CAS with extinction and ice water content derived from 2D- 
C and counterflow virtual impactor measurements. The con- 
clusion drawn from this analysis was that the former probes 
are measuring a combination of natural small crystals and 
shattered large ice crystals. Recently, Davis et al. (2009) 
compared optical depth derived from a number of instrument 
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probes (including CAPS) used during the WB57 Mid Lati- 
tude Cirrus Experiment (M1DCIX) with remote-sensing re- 
trievals from the Moderate Resolution Imaging Spectrome- 
ter (MODIS). In general the in situ probe optical depths were 
higher than the MODIS optical depths, and the CAPS opti- 
cal depths were the highest of the in situ probes. Davis et al. 
(2009) concluded that these results were consistent with past 
studies suggesting shattering artifacts in the CAS measure- 
ments. 

Despite evidence indicating that shattering at inlets may 
produce small-crystal artifacts, some recent publications 
have still argued for the common presence of high concentra- 
tions of small crystals in cirrus, and the measurements have 
been used at face value to make inferences about ice nucle- 
ation processes and cloud radiative properties. For example, 
Garrett et al. (2005) analyzed CAS measurements in sub- 
tropical anvil cirrus and concluded that small crystals dom- 
inate the crystal size distribution and anvil radiative prop- 
erties. Other studies have used the in situ measurements to 
suggest anomalous nucleation of crystals at cloud boundaries 
or evaporation zones (e.g. Cotton and Field, 2002). 

Karcher and Strom (2003) reported CVI measurements 
indicating cirrus clouds typically contained high concentra- 
tions (1-10 cm -3 ) of very small (diameter <20 /um) ice crys- 
tals. They showed that these large ice concentrations could 
be explained by ice nucleation driven by rapid gravity-wave 
driven cooling. However, the simulations only represented 
ice concentrations at the time and location where nucleation 
occurred. Cloud aging processes (entrainment, sedimenta- 
tion, and aggregation) decrease ice concentration with time. 
One might expect that the CVI measurements included aged 
cirrus or precipitation zones where the existence of numer- 
ous small crystals would not be expected. CVI measure- 
ments of ice concentration can be subject to shattering ar- 
tifacts if multiple shattering fragments contain sufficient sol- 
uble or insoluble material to produce detectable residual par- 
ticles. Hoyle et al. (2005) conducted a similar study com- 
paring parcel-model simulations with ice concentration mea- 
surements from a Multiangle Aerosol Spectrometer Probe 
(MASP, similar to CAS) and a CVI in cirrus during the Sub- 
sonic Aircraft: Contrail and Cloud Effects Special Study 
(SUCCESS). They also concluded that high ice concentra- 
tions (1-10 cm -3 ) are commonplace in midlatitude cirrus. 

Fridlind et al. (2004) used CAPS measurements as evi- 
dence for the importance of middle and upper tropospheric 
aerosols in the production of ice crystals that dominate the 
ice concentration in subtropical anvils. This issue is impor- 
tant for evaluating the relationship between aerosol proper- 
ties and cumulonimbus anvil microphysical properties and 
radiative effects. If boundary-layer aerosols dominate pro- 
duction of ice crystals in the clouds, then local pollution 
sources could affect anvil properties. If aerosols entrained 
into updrafts in the free troposphere dominate, then long- 
range transport and remote sources must be considered. As 
discussed further below, the conclusions of this study may 


need to be revisited given the indications of shattering arti- 
facts in CAPS measurements of small crystals in anvil cirrus. 

Recently, the Cloud Droplet Probe (CDP) was used to 
sample tropical anvil cirrus during the Tropical Warm Pool 
- Ice Crystal Experiment (TWP-ICE). Since it does not in- 
gest ice crystals into the sample volume through an inlet, 
but instead has probe arms, the CDP was expected to pro- 
duce relatively fewer shattering artifacts than the probes with 
inlet tubes. McFarquhar et al. (2007) compared CDP and 
CAS anvil size distributions and showed that the CAS con- 
centrations of small crystals appear to be overestimated by 
approximately two orders of magnitude when large crystals 
are present. The correlation between the ratio of CAS and 
CDP small-crystal concentrations and the concentration of 
large (>100/xm) crystals is consistent with the hypothesis 
that the enhanced CAS concentrations result from shattering 
of large crystals impacting the instrument inlet. 

The need to avoid the shattering problem and provide 
reliable information about small-crystal concentrations has 
motivated development of balloon-borne replicators. These 
replicators also provide sufficient resolution to determine ice 
crystal habits. Replicators were successfully launched dur- 
ing several field experiments (Miloshevich and Heymsfield, 
1997). In addition, the Video Ice Particle Sampler (VIPS) is 
designed to measure ice particles as small as ~ 10-20 /tm by 
impacting them on a transparent moving belt, followed by 
imaging the crystals with a video microscope (Heymsfield 
and Miloshevich, 1995). 

Here, we present CAPS measurements of ice crystal size 
distributions in tropical anvil cirrus from the the Tropical 
Composition, Cloud, and Climate Coupling (TC4) experi- 
ment. As in other datasets from instruments with inlets, the 
CAPS measurements indicate the presence of high concen- 
trations of small crystals throughout the vertical and horizon- 
tal extents and ages of the clouds sampled. We also compare 
the TC4 CAPS measurements with measurements from the 
2 Dimensional Stereo (2D-S) probe, which has the demon- 
strated ability to image ice crystals as small as 10 /xm (Law- 
son et al., 2006b). We find that the anvil cirrus small-particle 
concentrations indicated by the CAS component of CAPS 
generally exceed those indicated by 2D-S by 1-2 orders of 
magnitude in the 10-50 /mi size range where the instruments 
overlap. The CAS/2D-S concentration ratio is strongly corre- 
lated with large-crystal mass, suggesting that the discrepancy 
is caused by shattering of large crystals on the CAS inlet. 

Using the 2D-S dataset, we show that small crystals have 
little effect on extinction, radiative forcing, and radiative 
heating for the tropical anvil cirrus sampled during TC4. We 
show that high concentrations (as large as 30 cm -3 in tur- 
rets and >1 cm -3 in fresh convective anvil outflow) of small 
crystals do occur, but even in these locations, larger crystals 
dominate the cloud extinction. Also, it appears that the nu- 
merous small crystals detrained from deep convection do not 
persist as the anvil cirrus ages. Lastly, we show that con- 
centrations can exceed 1 cm -3 in isolated locations in aged 
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anvil cirrus well away from active convection; there is some 
evidence suggesting that the latter occurrences are associated 
with either small-scale convection or gravity waves. 

2 Instrument descriptions 

The TC4 mission, based in San Jose, Costa Rica, included 
extensive sampling of tropical anvil cirrus, using CAPS and 
2D-S measurements both on the WB-57 and the NASA DC- 
8 aircraft. CAPS consists of two instruments: the Cloud 
Aerosol Spectrometer (CAS) retrieves particle sizes by mea- 
suring light scattered off the particles and has an effective 
size range of 0.5-50 /rm; the Cloud Imaging Probe (CIP) is 
similar to 2D-C instruments used for decades and captures 
images of ice crystals with sizes ranging from about 25 to 
1550 /im. (The CIP pixel size is 25 /xm). CAS measures both 
forward and backward scattering. Only the forward scatter- 
ing signal is used in the data presented here, and in this mode, 
CAS is very similar to FSSP instruments commonly used for 
measurement of cloud and aerosol size distributions (Baum- 
gardner et al., 2001). 

In most previous field campaigns (including CRYSTAL- 
FACE and TWP-ICE), the CAS instrument has been flown 
with a shroud in front of the inlet. The purpose of the shroud 
is to straighten the airflow and remove dependence on air- 
craft angle of attack. However, the shroud provides an addi- 
tional surface for potential shattering of large crystals, and 
calculations show that small-crystal fragments coming off 
the shroud could get swept into the CAS sample volume (Mc- 
Farquhar et al., 2007). In fact, McFarquhar et al. (2007) ar- 
gued that the shroud was predominantly responsible for the 
shattering artifacts in the TWP-ICE CAS data. In the hope 
of at least partially mitigating the shattering problem, dur- 
ing TC4 the CAS probe on the DC-8 was flown without the 
shroud. 

The 2D-S probe is actually two independent, identical op- 
tical array instruments housed in one canister. Hydro meteors 
are imaged in the air flowing between the four detector arms 
(see Fig. 1). The 2D-S imaging instrument operates on prin- 
ciples similar to the CIP probe. However, it has significantly 
smaller pixel size (10 \x m) and faster electronics, resulting in 
size distributions extending down to 10 /xm (Lawson et al., 
2006b, 2008). It has been demonstrated that with the faster 
electronics the 2D-S can detect and count small particles at 
airspeeds up to ~230 ms -1 . Likewise, the smaller pixel size 
allows detection of crystals that might have been missed by 
the CIP and 2D-C probes with larger pixel size. The sizing of 
crystal images that are out of focus in both the CIP and 2D-S 
datasets are corrected based on diffraction theory (Korolev, 
2007). 

The 2D-S probe shares the problem of rapidly decreasing 
sample volume with decreasing particle size that CIP suffers 
from, and this decreasing sample volume is accounted for 
in the same way for each of the probes. The depth of field 



Fig. 1. 2D-S probe. The saucers on the detector arms are designed 
to minimize shattering fragments from entering the sample volume, 
which is near the middle of the path between the detector arms 
for small crystals. However, shattering fragments do still reach the 
sample volume and are apparent as clusters of small particles with 
very small inter-arrival times. These artifacts are removed using a 
detailed statistical analysis, as described by Baker et al. (2009) for 
the case of splashing raindrops. See text for details. 


is assumed to be proportional to the square of the particle 
size. Experiments with glass beads were used to determine 
the coefficient in this relationship. A thorough investigation 
of the errors associated with depth of field effects on sample 
volume and corrections for out of focus particles has not been 
done. Therefore it is difficult to quantify the uncertainties in 
2D-S measurements of small-crystal concentrations. 

Although shattering artifacts do still occur in instruments 
such as CIP and 2D-S as a result of large crystals impacting 
the detector arms (Field et al., 2006), the shattering problem 
is expected to be less severe than for the CAS probe with 
its inlet. One reason for this expectation is that the detec- 
tor arms protruding into the airstream are further from the 
sample volume than is the case for the CAS inlet. In addi- 
tion, the 2D-S probe has sharp disks on the interior of the 
detector arms that are intended to deflect shattering artifacts 
away from the sample volume. While fewer shattering ef- 
fects are expected with the construction of the 2D-S com- 
pared to the CAS or FSSP, numerous shattering fragments 
do still reach the sample volume and are apparent as groups 
of closely spaced particles. These spurious effects are re- 
moved (assuming ice particles are randomly distributed in 
clouds) using algorithms described by Baker et al. (2009) for 
the case of splashing raindrops, with the exception that tests 
for image roundness are excluded when processing ice data. 
The main strategy uses the close spacing of spurious particles 
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to eliminate them, but some image processing techniques are 
also used. If the filtering algorithms are not employed the 
shattering effects on the 2D-S can be great at the smaller 
sizes when the natural concentration of small particles is low 
(Baker et al., 2009). The DC-8 CIP data were processed with 
a similar algorithm for removing shattering artifacts based on 
particle interarrival times. 

It should be noted that if a large crystal shatters on the 
probe housing and only one of the fragments passes through 
the sample volume, then the interarrival time analysis will 
not identify the sampled fragment as an artifact. Although 
we believe such events are not common, they could result 
in an artificial enhancement of the 2D-S small-crystal mea- 
surements. Another possible source of artifacts is small crys- 
tals that bounce off the probe arms and are deflected into the 
sample volume. If either of these processes are significant, 
then the natural small crystal concentration will actually be 
less than the 2D-S results presented below. Our primary con- 
clusion - that small crystal contribute relatively little to the 
TC4 anvil cirrus radiative properties - would thus not be af- 
fected by these artifacts. Also, if these artifacts are signifi- 
cant in the 2D-S measurements, it would imply that the CAS 
measurements of small -crystal concentration are exaggerated 
even more than we suggest below. 

On the WB-57, the 2D-S probe was mounted under the 
wing about 2 m back from the wing leading edge, with the de- 
tector axis about 16 cm below the skin of the wing. There was 
concern that shattering artifacts resulting from large crystals 
impacting the wing leading edge might contaminate the in- 
strument in this mounting location. This issue will be dis- 
cussed further in Sect. 2.4. On the DC-8, the 2D-S and CAPS 
probes were mounted on wing pylons a couple of feet be- 
low the wing with the detectors slightly forward of the wing 
leading edge, and shattering from airframe surfaces should 
not have been a problem. The results presented here focus 
primarily on the DC-8 measurements. 

2.1 TC4 anvil cirrus measurements 

During the July-August 2007 TC4 mission, tropical anvil cir- 
rus with a range of ages and generated by a variety of types 
of deep convection were sampled, including fresh and aged 
outflow from Inter-Tropical Convergence Zone (ITCZ) sys- 
tems, aged anvil cirrus generated by strong systems in the 
Panama Bight, and fresh outflow from systems along the 
coast of Costa Rica. In addition, on two DC-8 flights, the 
aircraft flew through the tops of active turrets in unusually 
intense ITCZ systems. It should be noted that most of the 
anvil cirrus sampling was conducted well away from active 
convection for safety reasons. 

First, we address the issue of shattering artifacts by com- 
paring the CAPS measurements with the 2D-S measurements 
in the overlapping size range. As discussed above, 2D-S 
is expected to be less susceptible to shattering, and shatter- 
ing artifacts are removed with the post-processing interarrival 


time analysis. As also discussed above, the shroud in front 
of the CAS probe on the DC-8 was removed for the TC4 
campaign, potentially substantially reducing the frequency 
of shattering artifacts. Hence, it was anticipated that the con- 
trast between the small-crystal concentrations indicated by 
the CAS probe with its inlet and the 2D-S instrument might 
be reduced compared to the TWP-ICE results presented by 
McFarquhar et al. (2007), where the CAS concentration typ- 
ically exceeded CDP by approximately two orders of magni- 
tude. 

On the flights of 8 August, the aircraft sampled anvil cirrus 
over the ocean south of Costa Rica, and the DC-8 and WB- 
57 pilots were instructed to successively fly through the same 
cloud at the same altitudes. Figure 2 shows a comparison 
between the CAPS (CAS and CIP) and 2D-S measurements 
from DC-8 and WB-57 flight legs at ~ 1 2 k m and ~ 11.4km. 
The agreement between size distributions derived from the 
CIP and 2D-S imaging instruments is excellent throughout 
their overlapping size range (>50 /zm). The 2D-S measure- 
ments from the two aircraft agree reasonably well, particu- 
larly at sizes larger than ~150/zm. The WB-57 2D-S mea- 
surements indicate about 2-3 times as many smaller crystals 
than indicated by the DC-8 2D-S. This difference may be 
an indication of shattering off the leading edge of the WB- 
57 wing affecting the 2D-S measurement at the wing-hatch 
location. Such artifacts would be spatially separated in the 
sheared flow between the wing leading edge and the instru- 
ment and would not be identified in the interarrival time anal- 
ysis. The issue of shattering on the WB-57 wing affecting 
the 2D-S measurements under the wing is addressed further 
in the next section. 

In contrast to the agreement between the imaging instru- 
ments, the CAS measurements indicate about an order of 
magnitude higher concentrations of small ( <50 //m) crys- 
tals than do the 2D-S measurements on the 12-km DC-8 
flight leg. At 11.4km, where more large crystals are present, 
the CAS concentrations are about two orders of magnitude 
higher than the 2D-S concentrations. Although there is more 
uncertainty in the 2D-S concentrations in the 10-30 fi m size 
range than at large sizes (due to depth of field uncertain- 
ties, statistical uncertainties resulting from small sample vol- 
umes, and uncertainties due to imperfect shattering removal), 
the CAS/2D-S discrepancy extends out to the maximum size 
limit of the CAS probe (~50jU.ni). Comparisons between 
CAS and 2D-S on other flights indicate consistent results: 
when the natural concentration of small crystals is relatively 
low (i.e., a few hundreds per liter or less) as is the case most 
of the time in the anvil cirrus, and when the concentration of 
large ice is relatively high (order 0.001 cm -3 of 500 /zm crys- 
tals), the ratio of CAS/2D-S concentration measurements is 
typically 10 to 100. 

For comparison with the anvil cirrus measurements, we 
also show in Fig. 2 2D-S and CAPS measurements in driz- 
zling stratocumulus sampled during a boundary-layer run 
at about 1 km. In this case, there are relatively few large 
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TC-4 8 Aug DC-8 anvil cirrus at 12 km 



TC-4 8 Aug DC-8 anvil cirrus at 1 1 .4 km 



TC-4 6 Aug DC-8 drizzle at 1 km 



Fig. 2. Size distributions averaged over flight legs when WB-57 and 
DC-8 successively flew along the same flight legs at approximately 
12 km (top) and 11.4 km (bottom) on 8 August are shown in the 
top two panels. The DC-8 CIP distribution agrees well with the 
2D-S dish'ibution; however, the concenhations of small particles 
indicated by CAS are 1-2 orders of magnitude higher than indicated 
by 2D-S. For comparison with the anvil measurements, we show 
2D-S and CAPS size distributions from measurements in drizzling 
stratotocumulus where shattering (or splashing) should have little 
impact and the small-droplet concentrations measured between the 
two probes are in reasonable agreement. 

droplets and splashing should have a minimal impact on the 
measured concentrations of small droplets. In fact the inter- 
arrival time analysis used in the 2D-S data processing indi- 


cates very few splashing artifacts. Although there are some 
differences in sizing between CAS and 2D-S in the over- 
lapping size range (<50 /xm), the concentrations are similar 
(particularly in contrast to the comparisons in anvil cirrus). 

It is worth noting that in areas with relatively high natu- 
ral small ice concentrations, effects of shattering on the CAS 
do not have as large a relative impact on ice concentrations 
as they do in regions with low natural ice concentrations. 
A simple example will explain why this result is expected: 
assume that a fixed mass of crystals large enough to gener- 
ate shattering artifacts produces 1 cm~ 3 enhancement in the 
measured small-crystal concentration. If the natural concen- 
tration of small crystals were 0.1 cm -3 , then the relative en- 
hancement caused by shattering would be large (about an or- 
der of magnitude); whereas if the natural concentration of 
small crystals were greater than 1 cm -3 , then the relative im- 
pact of shattering on measured ice concentration would be 
less important (about a factor of 2). This point is exempli- 
fied by comparison of the particle size distributions in Fig. 3, 
which shows CAPS and 2D-S measurements on 22 July in 
regions with high and low natural small ice abundance , in 
conjunction with approximately the same amount of large 
ice. Note, however, that the large ice only extends out to 
500 /xm, compared with 1 to 3 mm in Fig. 2. In Fig. 3, the 
2D-S size distribution is shown with and without removal of 
shattering artifacts identified with the interarrival time anal- 
ysis. In the presence of low natural ice, shattering increases 
the small ice concentration detected by the 2D-S probe by 
about an order of magnitude, and the CAS concentration is 
enhanced by as much as 2 orders of magnitude. In contrast, 
for the case with numerous natural small particles, shatter- 
ing enhances the 2D-S small-crystal concentration by only 
a factor of 2, and the CAS measurements agree reasonably 
well with 2D-S. Although the number of shattered particles 
is on the order of 0.1 cm -3 in each case, the relative effect 
is much greater in the case with few natural small crystals. 
Note that the two regions are only 10 km apart near the top 
of aged anvil cirrus. 

If the discrepancy between the CAS and 2D-S concentra- 
tions is caused by shattering on the CAS inlet, one would 
expect a correlation between the CAS/2D-S concentration ra- 
tio and the amount of large-ice mass available for shattering. 
Figure 4 shows that such a correlation does exist. For the 
three hours of anvil sampling on 8 August, the CAS/2D-S 
ratio and ice mass in crystals larger than 175 /xm are strongly 
correlated (Pearson correlation coefficient of 0.74). The ice 
water content is calculated from the 2D-S images using the 
area-to-mass parameterization described by Baker (2006). 
Our conclusion from this analysis is that even with the shroud 
removed, when the natural concentration of small crystals 
is on the order of 0.1cm -3 or less and 500 /xm ice is at 
least about 0.001 cm -3 , the small-particle concentrations in- 
dicated by the CAS probe appear to still be dominated by 
shattering artifacts. 
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Fig. 3. CAPS and 2D-S size distributions (with and without shattering artifacts removed) are shown for two cases. Left panel: a case with 
numerous natural small crystals for which shattering has a relatively small impact on small-crystal concenhation. Right panel: a case with 
low natural small-crystal concentrations for which shattering artifacts significantly affect the measurements. 


2.2 Small crystal concentration - large crystal mass cor- 
relations 

Another way to investigate the impact of ice crystal shatter- 
ing artifacts is to examine the correlation between small- and 
large-crystal abundance. In convective updrafts or very fresh 
outflow, the concentration of both large and small crystals 
may be predominantly controlled by dilution of the cloud 
as ambient air is entrained, in which case the concentra- 
tions of small and large crystals will likely be correlated. On 
the other hand, as anvil cirrus ages, differential gravitational 
settling should separate the large and small crystals. One 
would expect to find places in the cloud where large crystals 
are present without small crystals and vice versa. For ex- 
ample, at the base of the cloud, the size distribution should 
be dominated by large crystals that sedimented from above, 
with relatively few small crystals present. Breakup of sub- 
limating large ice crystals could produce some smaller frag- 
ments (Bacon et al., 1998), however this process would not 
likely produce large concentrations of small crystals. If the 
measured small-crystal concentration is dominated by shat- 
tering artifacts, we would expect to see a strong correlation 
between small and large crystal abundance in the measure- 
ment whenever large crystals are present, even in aged anvil 
cirrus, and even in precipitation zones where natural small 
crystals are actually absent. Assuming that the sizes of frag- 
ments generated by large-crystal shattering are roughly inde- 
pendent of the large-crystal size, the number of small frag- 
ments generated should be proportional to the total mass of 
large crystals that shatter. Such correlations were found in 
the CRYSTAL-FACE and other datasets (Heymsfield, 2007; 
McFarquhar et al., 2007). 


8 August anvil cirrus 
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Fig. 4. Scatterplot of the DC-8 CAS/2D-S concentration ratio ver- 
sus ice water content in large (>175^111) ice crystals in the anvil 
cirrus sampled on 8 August. The strong correlation suggests that 
the CAS concentrations are exaggerated as a result of shattering ar- 
tifacts. 

Scatterplots of measured small-crystal concentration ver- 
sus large-crystal mass for the 8 August anvil cirrus are plot- 
ted in Fig. 5. Again, we restrict ourselves to the flight legs 
at 1 1.4 and 12 km such that the DC-8 and WB-57 were sam- 
pling aged anvils with comparable crystal size distributions. 
For these plots, we have used the same size ranges for small 
(10-55 /zm) and large (> 1 20 /cm) crystals for the different 
instruments. In the DC-8 CAPS dataset (Fig. 5a), small- 
crystal concentration and large-crystal mass are strongly 
correlated throughout the aged anvil cirrus. Note also the 
very high concentrations of small crystals in the CAS data. 
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Fig. 5. Scatterplots of small-crystal concentration versus large- 
crystal mass are plotted for the 8 August anvil measurements. We 
restrict this analysis to the time periods when the DC-8 and WB-57 
were flying in the 1 1.4-12 km altitude range. The three panels cor- 
respond to measurements from the DC-8 CAPS, DC-8 2D-S, and 
WB-57 2D-S. Similar correlations are apparent in anvil cirrus sam- 
pled on other days. See text for discussion. 


DC-8 2D-S 8 Aug anvil cirrus 



DC-8 CAPS 8 Aug anvil cirrus 
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with values often exceeding 1 cm' 3 and occasionally ex- 
ceeding 10 cm' 3 . Similar strong correlations between small- 
crystal concentration and large-crystal mass (not shown) are 
apparent in the CRYSTAL-FACE and TWP-ICE CAPS anvil 
cirrus datasets, again suggesting that, in the presence of large 


ice when the natural concentration of small crystals is rela- 
tively low, the small-crystal concentration measured by CAS 
is dominated by shattering artifacts. 

In contrast, the DC-8 2D-S scatterplot (Fig. 5b) is more in 
line with what we would expect for aged anvil cirrus from 
physical considerations. There is a relatively weak corre- 
lation between small and large crystal abundance, and the 
dataset does include locations with plentiful large crystals 
but very few small crystals. Note also that the DC-8 2D-S 
small-crystal concentrations are much lower than the CAPS 
measurements suggest, with values only very rarely exceed- 
ing 1 cm' 3 . 

Looking next at the WB-57 2D-S scatterplot (Fig. 5c), we 
see a stronger correlation between small-crystal concentra- 
tion and large-crystal mass than apparent in the DC-8 2D-S 
dataset, as well as an absence of occurrences with consid- 
erable abundance of large crystals without numerous small 
crystals. As discussed above, it was recognized in the plan- 
ning stages of TC4 that the mounting location of the 2D-S 
probe on the WB-57 was less than ideal. With the probe on a 
wing hatch toward the rear of the wing, there was a possibil- 
ity of contamination from small-crystal fragments generated 
by large crystals impacting the wing leading edge and shat- 
tering. Fragments from wing leading edge shattering events 
cannot be readily removed by examining particle interarrival 
times, because the fragments get spatially separated in the 
sheared flow between the wing leading edge and the probe lo- 
cation 2 m back on the wing. This analysis suggests that such 
contamination may indeed have affected the WB-57 2D-S 
size distribution measurements. To be sure the 2D-S results 
presented are not contaminated by shattering artifacts, in the 
remainder of this paper we will, unless otherwise noted, fo- 
cus on the DC-8 2D-S measurements. We emphasize that the 
problem of apparent shattering artifacts in the WB-57 2D- 
S data is an airframe and mounting location problem rather 
than an instrument problem. 

2.3 Quantifying shattering artifacts 

The comparisons shown above suggest that shattering of 
large crystals can be a large source of artifacts in the CAS 
measurements. Field et al. (2003) and Field et al. (2006) 
provided a quantitative measures of the impact of shatter- 
ing on FSSP and imaging probes (e.g., 2D-C) measurements, 
respectively, by reporting the ratio of measured ice concen- 
trations to those with shattering artifacts removed (using in- 
terarrival times). However, this approach does not allow 
users of other FSSP datasets to estimate the concentration of 
shattering artifacts under different conditions. As discussed 
above, the relative impact of shattering will depend both on 
the number artifacts and the concentration of natural small 
crystals. For example, Arctic mixed-phase clouds have much 
lower natural ice concentrations (typically a few per liter; 
Fridlind et al., 2007) than anvil cirrus (typically on the or- 
der of 0.1 cm -3 ). Thus, the relative impact of shattering on 
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ice concentrations measured by 2D-S probes may be much 
more severe for the Arctic clouds than indicated by Field 
et al. (2006) for anvil cirrus. 

As a means of quantifying the effect of shattering on the 
CAS measurements, we present here the correlation between 
spurious CAS concentration and the 2D-S large-crystal mass. 
We have selected averaging periods with reasonably homo- 
geneous cloud conditions by applying the following criteria: 
(1) continuous cloud sample of at least 10 s, (2) 2D-S con- 
centration >50 /xm to remain within 0.5 and 1.5 times the 
average concentration over the averaging period, (3) 2D-S 
concentration <50 /xm to remain within 0.5 and 1.5 times 
the average concentration over the averaging period, and 
(4) average concentrations of 2D-S particles >50 /xm to be 
>10 -6 cm -3 . We also exclude time periods for which the 
CAS concentration is less than 10 times greater than the 
2D-S concentration over their overlapping size range (~12- 
50 /xm) such that the discrepancy between CAS and 2D-S 
can be assumed to be shattering artifacts. Using these re- 
stricted datasets for different flights, we compute linear least- 
squares regression fits. 

The resulting scatterplots of spurious CAS (> 12 /xm) con- 
centration, extinction, and IWC versus 2D-S IWC and the 
linear fits are shown in Fig. 6 for anvil cirrus sampled on the 
flight of 5 August. With the dataset restricted to the relatively 
homogeneous time segments, the correlations are excellent 
(e.g., a Pearson correlation coefficients >0.97). It should be 
noted that the slope of the correlations varies significantly 
from flight to flight, and the fits presented in Fig. 6 should 
be taken as a very rough guide to the quantitative impact of 
shattering on the CAS probe. 


3 Implications of evidence for shattering artifacts 

The evidence presented here suggests that small-crystal con- 
centrations measured by the CAS probe (with or without 
the shroud in front of the inlet) are dominated by shattering 
artifacts when large crystals are present and natural small- 
crystal concentrations are relatively low. The great majority 
of the anvil cirrus measurements made during CRYSTAL- 
FACE had considerable amounts of large ice, and thus, when 
the natural small ice concentration is relatively low, the CAS 
concentrations were probably dominated by shattering arti- 
facts, as they appear to be in the TWP-ICE and TC4 anvil cir- 
rus datasets. Garrett et al. (2005) analyzed CRYSTAL-FACE 
CAPS data and concluded that in anvil cirrus sampled on 21 
July 2002, crystals smaller than 50 /xm dominated the size 
distributions and radiative properties. Baumgardner et al. 
(2005) used CRYSTAL-FACE CAS data to infer habits of 
small crystals, and Chepfer et al. (2005) compared remote- 
sensing retrievals of ice crystal habits in CRYSTAL-FACE 
anvil cirrus with habit information retrieved from CAS mea- 
surements. We suggest that the conclusions of these studies 
should be revisited given the likelihood that CAS was often 
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Fig. 6. Scatterplots of spurious CAS concentration (top), extinction 
(middle ), and and IWC (bottom) (1 1-50 /xm) versus 2D-S ice water 
content for time segments with relatively homogeneous clouds (see 
text). The lines shown are linear least-squares regression fits, and 
the slopes are given in the legends. We caution that these correla- 
tions should be taken as very approximate estimates to the magni- 
tude of shattering artifacts that might be seen in other clouds sam- 
pled with CAS. 


primarily measuring shattering artifacts in the anvil cirrus. 
The WB57 Mid Latitude Cirrus Experiment (MIDCIX) in- 
cluded sampling of a wide range of midlatitude cirrus types, 
some of which may have been sufficiently lacking in large 
crystals such that the CAS data may be valid, and some of 
which may have been contaminated by shattering. In any 
case, we suggest these datasets should be used with caution. 
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The indications of shattering contamination in the TC4 
CAS measurements, with the shroud removed, may have 
even broader implications. The CAS probe without the 
shroud has an inlet geometry very similar to the FSSP in- 
struments used for the past two decades to measure ice crys- 
tal size distributions in cirrus, including tropical anvil cirrus. 
However, direct comparisons between FSSP and instruments 
such as 2D-S and CDP, along with use of FSSPs with inter- 
arrival times recorded, are required to further evaluate un- 
der what conditions the historical dataset can be relied upon. 
Nevertheless, the evidence presented here suggests that, in 
general, it would seem prudent to be suspicious of ice crystal 
concentrations inferred from any probes with inlets or sur- 
faces protruding into the airstream (without the use of in- 
terarrival times to identify and remove shattering artifacts), 
including CVIs, whenever large crystals are present. 

Fridlind et al. (2004) used CRYSTAF-FACE CAPS mea- 
surements in subtropical anvil cirrus as evidence for the 
predominance of free-tropospheric aerosols entrained into 
the deep convective updrafts (rather than boundary-layer 
aerosols in the inflow region) for production of ice crystals 
that ultimately get detrained into the anvils. Cloud-resolving 
model simulations were run with all aerosols in the verti- 
cal column, with aerosols above 6 km removed, and with 
aerosols above 2 km removed. Only when all aerosols in 
the vertical column were included did the simulations re- 
produce the size distributions measured by CAPS in the up- 
per parts of the anvils, with a distinct mode in the num- 
ber distribution at around 20-30 // m and very large concen- 
trations of small crystals. If aerosols above 2 km were re- 
moved in the simulation, the mode in the number distribution 
shifted out to ~60 /xm, and the small crystals were largely 
absent. However, given the results presented by McFarquhar 
et al. (2007), the results presented here, and the fact that, in 
general, there were considerably more large crystals in the 
CRYSTAF-FACE anvil cirrus than in the cirrus sampled dur- 
ing TWP-ICE and TC4, it seems very likely that the small- 
crystal mode in the CAPS data was strongly influenced by 
shattering artifacts. 

Although other lines of evidence for the predominance of 
mid-tropospheric aerosols were presented by Fridlind et al. 
(2004) involving measurements made in a mixed-phase re- 
gion of a convective updraft, the key issue of interest for the 
impact of cumulonimbus clouds on climate is what controls 
the size distribution of ice crystals detrained into the anvil 
cirrus. If the CAPS small-crystal concentrations are indeed 
dominated by shattering artifacts, then one must conclude 
that the agreement between the simulated ice crystal size 
distribution with the CAPS distributions indicates the simu- 
lated anvil had excessive numbers of small crystals. There 
are a number of physical processes in the model, such as 
entrainment and ice crystal aggregation, that are not well 
constrained, and misrepresentations of these processes might 
well result in excessive ice concentrations. The spatial res- 
olution (500 m) used in the simulations is inadequate for re- 


solving the eddies that dominate entrainment into convective 
updrafts. An experiment was done with 250 m spatial reso- 
lution, but even finer resolution may be required. 

Other evidence from the CRYSTAF-FACE mission points 
to the possibility that anvil ice crystals predominantly orig- 
inated from aerosols in the lowermost troposphere. Cziczo 
et al. (2004) analyzed the chemical composition of residual 
aerosols from anvil ice crystals. They found that 70% of the 
residual particles were sea salt, mineral dust, or fly ash. None 
of these particle types were present in significant concentra- 
tions in clear free -tropospheric air, suggesting the source was 
the planetary boundary layer. It should be noted that this 
analysis only included residual particles larger than about 
0.2 fim diameter, and smaller aerosols probably contribute 
significantly to droplet activation in strong cumulonimbus 
updrafts. Twohy and Poellot (2005) analyzed residual par- 
ticles larger than 0. 1 /xm collected with a counterflow virtual 
impactor. They found that about 37% of the small residual 
nuclei and 49% of the large residual nuclei were either sea 
salt or industrial particles, again suggesting a strong depen- 
dence of anvil ice concentrations on boundary-layer aerosols. 
This discussion underscores the fact that accurate measure- 
ments of ice crystal size distribution are required for evalu- 
ation of aerosol impacts on cirrus microphysical properties 
and radiative effects. 

4 Contribution of small crystals to tropical anvil cirrus 
radiative properties 

As discussed above, past measurements using instruments 
with inlets have generally indicated that small crystals often 
dominate, or at least contribute significantly to, the surface 
area density (and accordingly, extinction) in cirrus clouds. 
Here, we present several cases from the TC4 mission using 
2D-S data that suggest, at least for the anvil cirrus sampled 
during this mission, small crystals contribute relatively little 
to the extinction. 

An unprecedented level of aircraft coordination was 
achieved during TC4, with the ER-2 and either one or both 
of the WB-57 and DC-8 aircraft flying in a stacked forma- 
tion for many hours on several flights. During much of the 
coordinated-flight time, the WB-57 and DC-8 were flying in 
the anvil cirrus. Besides providing a great deal of in situ data 
for comparison with cloud remote-sensing retrievals from 
various instruments on the ER-2, the aircraft coordination 
also allows us to use the remote-sensing information to pro- 
vide context for the in situ measurements. In particular, we 
can tell where the aircraft measurements were made relative 
to cloud top and cloud base. 

For cases with significant periods of good aircraft coor- 
dination, we present the merged ER-2 Cloud Physics Fidar 
(CPF) and Cloud Radar System (CRS) product time cur- 
tain (i.e., retrieved cloud property versus time and height 
along the ER-2 flight track), and we indicate on these curtain 


Atmos. Chem. Phys., 9, 5519-5537, 2009 


www.atmos-chem-phys. net/9/55 19/2009/ 



E. J. Jensen et al.: Small Crystals in Anvil Cirrus 


5529 


Merged CPL/CRS cloud product 8 Aug ITCZ anvils 
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Fig. 7. Top panel: curtain plot the backscatter derived from a combination of the ER-2 CPL and CRS measurements versus time and height. 
The black trace indicates the altitude of the DC-8 flying below the ER-2; where the altitude trace is thick, the horizontal separation between 
the DC-8 and ER-2 was less than 5 km. Middle plot: time series of surface area density, ice concentration, ice water content and effective 
radius, all derived from the 2D-S measurements. Bottom plots: number and surface area distributions averaged over the four constant- 
altitude flight legs indicated in the top panel. Leg averaged ice concentrations, effective radii, small-crystal (10-50 /im) surface area density 
and large-crystal (>120 gm) SAD (fimr cm -3 ) are indicated in the legends. 


plots the altitude of the DC-8 flying below the ER-2. The 
merge product contains coincident vertical curtains of CPL 
backscatter and CRS reflectivity, which when combined, 
gives the full vertical cloud distribution not attainable with 
either instrument alone. Merged signal images contain nor- 
malized values (0-1) of lidar backscatter ratio using the log 
of the ratios from 0.5 to 500.0 and normalized values (0-1) 
of radar reflectivity (dBZ) using the range —30 to 20dBZ. 
The type of signal making up the cloud parcels (i.e., CPL, 
CRS, or both) is also tracked. Next, we present time series 
of bulk cloud properties calculated from the 2D-S measure- 
ments. Lastly, for time periods indicated on the CPL/CRS 
curtain plots, we present averaged 2D-S number and surface 
area density distributions. The average contributions to total 
projected surface area density (SAD) from small (10-50 gm) 


and large (>120 gm) crystals, as well as the average ice con- 
centrations and effective radii, are given in legends on these 
plots. We use the standard definition of effective radius here: 


^eff = 


3IWC 

2/Jext Pi 


(1) 


where IWC is the ice water content, /3 ex t is the extinction, 
and pi is the bulk ice density (0.9 g cm -3 ). 

We begin with the anvil cirrus sampled over the eastern 
Pacific, south of Costa Rica by the DC-8 and ER-2 on 8 Au- 
gust (Lig. 7). On this day, we achieved the greatest degree 
of aircraft coordination, with the aircraft vertically stacked 
(to within 5 km horizontal separation) for long time periods. 
As indicated by the altitude trace on the CPL/CRS curtain 
plot, for most of the anvil-sampling time period, the DC-8 
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Fig. 8. As in Fig. 7, but for anvil cirrus generated by relatively weak ITCZ convection sampled on 22 July. Since the flight pattern primarily 
consisted of ascents and descents through the clouds, we present size distributions averaged over altitude ranges rather than time periods. 


was well below the cloud top altitude of 13-13.5 km. Also, 
the clouds generally extended at least a few km below the 
DC-8 flight altitude. Only on the last leg at ~12 km was the 
DC-8 near cloud top. The time series plot shows the large de- 
gree of variability typically observed in cirrus. The anvil cir- 
rus sampled on this flight was generally relatively aged, with 
ice water contents less than 0.2 gm“ 3 and total ice concen- 
trations typically on the order of 0.1 cm -3 and almost never 
exceeding 1 cm -3 . The size distributions averaged over the 
four constant-altitude flight segments indicate clear evidence 
of gravitational size sorting in the anvils, with fewer large 
crystals and more small crystals as one ascends higher in the 
anvils. Ice crystals smaller than 55 /zm contribute very little 
(less than 6%) to the total surface area density at any vertical 
level in the anvils. 

The effective radii derived from 2D-S ice water content 
and extinction range from about 30 to 60 /zm deep in the 
cloud, but are limited to about 30-35 /zm near cloud top. 


These cloud-top effective radii are in good agreement with a 
variety of remote-sensing retrievals, including MAS, MODIS 
(Hong et al., 2007), and the Solar Spectral Flux Spectrometer 
(SSFR) (P. Pilewskie, personal communication, 2008). Note 
that we are not arguing that this comparison with remote- 
sensing measurements serves as strong evidence for the accu- 
racy of the 2D-S measurements. The uncertainties in remote- 
sensing retrievals of cloud properties and considerable and 
difficult to quantify, precluding their use as a baseline for 
evaluating other measurements. We simply note the consis- 
tency between the different measurements. 

We next present the measurements made in outflow anvils 
from typical eastern-Pacific ITCZ convection on 22 July 
(Fig. 8). In this case, the cloud tops were near the DC-8 
maximum altitude, and the aircraft sampled from cloud top 
(about 11.5 km) down to about 8 km in a “saw-tooth” pat- 
tern. The height dependence of effective radius is evident; 
for example as the aircraft descends during the 14.3-14.4 UT 
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24 July convective turret (11-12 km) 





Fig. 9. Time series of ice concentration, and IWC (top panel) and vertical wind speed (middle panel) are shown for an active turret sampled 
on 24 July. Number and surface area distributions are shown in the bottom panels. 


time period from about 12 to 8 km, R e g- increases from about 
3CM-0 to nearly 80 /xm. The size distributions (in this case 
averaged over height ranges) again show the vertical size 
sorting and the limited contribution of small crystals to SAD. 

On the DC-8 flights of 22 and 24 July, very fresh outflow 
from anomalously strong ITCZ convection was sampled, as 
well as the tops of active turrets at 1 1-12 km. Figure 9 shows 
time series and size distributions in the strongest convection 
sampled on 24 July. Remote-sensing data from the ER-2 is 
not shown for this case because we did not have close air- 
craft coordination; however, the nadir viewing Second Gen- 
eration Precipitation Radar (APR-2) (Sadowy et al., 2003) on 
the DC-8 provides information about the convective system 
structure and dynamics below the aircraft. Imagery from the 
APR-2 on board the DC-8 is shown in Fig. 10. The DC-8 first 
flew over a portion of aged outflow followed by a very exten- 
sive region of active convection including four well defined 
convective cells (U1-U4) and a wide transitional area (T). Of 
the four active cells, U2 was the most developed vertically, 
with the peak of vertical wind estimated at around 17 ms -1 
at 9 km (the top of the visible range for the radar), large par- 
tially frozen hydrometeors lofted at least to 8 km, and a cou- 
plet downdraft at lower levels, suggesting that this cell was in 


the late stages of maturity. The last cell (U4) shows a region 
of maximum updraft at about 8.5 km with numerous large 
hydrometeors peaking at 5.5 km, the absence of a downdraft 
and relatively depleted volume underneath suggest that this 
cell was at a younger stage than U2. 

Ice concentrations (Fig. 9) in this cloud were extremely 
high (average values of 5.3 cm -3 and a peak value of 
32 cm -3 ). Even with these very large numbers of small crys- 
tals present, the surface area density is still dominated by 
large crystals, with crystals smaller than 55 /xm contributing 
less than 5% to the total SAD. It is also interesting that the 
mean effective radii in these events (~47 /xm) are somewhat 
larger than found in the tops of mature anvil cirrus (~30- 
40 /xm). This observation attests to the ability of the strong 
updrafts (20 ms -1 at ~12km) to loft large ice crystals (grau- 
pel) into the upper troposphere. 

The measurements presented here suggest that small- 
crystal concentrations decrease dramatically, from up to 
30 cm -3 in turrets, to a few cm -3 in fresh outflow, to 
<0. 1 cm -3 in anvil cirrus with ages ranging from about 
30 min to a few hours. It might seem counterintuitive that as 
anvil cirrus ages and larger crystals sediment, the upper parts 
of the anvil don’t appear to become increasingly dominated 
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Fig. 10. Imagery from the APR-2 radar showing reflectivity (top) and Doppler velocity (bottom). See text for discussion. 


by small crystals. However, other physical processes, such 
as sublimation, deposition growth, and aggregation, tend to 
decrease the small-crystal concentration as the cloud ages. 
In particular, in moderately subsaturated (supersaturated) air 
crystals will sublimate (grow) rapidly through the small- 
crystal size range (Fig. 11). Ice saturation ratio measure- 
ments in the TC-4 anvil cirrus typically indicate a great deal 
of variability with values ranging between about 0.8 to 1.2. 
Under such conditions, small crystals will either sublimate 
completely or grow out of the small-crystal size range on 
time scales of about 10-60 min. The measurements indicat- 
ing a lack of small crystals suggest that these sinks for small 
crystals likely dominate over any potential sources of small 
crystals (e.g., nucleation of new crystals in supersaturated re- 
gions) in the aging anvils. 

The calculations presented in Fig. 1 1 assume a range of 
values for the deposition coefficient. If the value is much 
smaller than 0.1, as suggested by a recent laboratory mea- 
surement (Magee et al., 2006) (a^ep ^ 0.006), then the small 
crystals could persist for hours even with moderate subsat- 
uration or supersaturation (Fig. 11). Past studies comparing 
numerical simulations and measurements of ice concentra- 
tions in wave clouds and synoptic cirrus have suggested that 
the deposition coefficient cannot be much less 0.1 (Jensen 
et al., 1998; Comstock et al., 2008). The fact that the small 
crystals do not persist in the anvil cirrus is also inconsis- 
tent with deposition coefficients smaller than 0.1, which 
would prevent 10-50 /zm crystals from growing or sublimat- 
ing rapidly. 

Since small crystals contribute very little to extinction in 
the anvils, it seems unlikely that they affect other radiative 
properties, such as radiative forcing and radiative heating. In 



Fig. 11. The times required for an ice crystal to either completely 
sublimate (in subsaturated air) or double it’s size (in supersaturated 
air) plotted versus ice saturation ratio for three different assumed 
deposition coefficients. These calculations assumed a temperature 
of 220 K and a pressure of 200 hPa. Spherical particle growth is as- 
sumed. Given the relative humidity variability typically observed in 
anvil cirrus, the fact that the plentiful small crystals detrained from 
convection do not persist is inconsistent with deposition coefficients 
much less than ~0.1. 

order to confirm this hypothesis, we have run radiative trans- 
fer calculations with cloud properties specified based on the 
2D-S measurements. We use the 8 August DC-8 and WB-57 
2D-S measurements to construct ice crystal size distributions 
versus height (Fig. 12). As discussed above, there is evi- 
dence that the WB-57 2D-S concentrations of small crystals 
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are exaggerated because of shattering on the wing leading 
edge. As a result, we may be overestimating the effects of 
small crystals in our calculations. We use a two-stream ra- 
diative transfer model (Toon et al., 1989) with the temper- 
ature and water vapor profiles specified based on the DC-8 
dropsonde at 12.89 UT. Mie scattering is used to calculate 
the cloud single-scattering properties. 

Figure 13 shows the top-of-atmosphere solar, infrared, and 
net radiative forcing versus the minimum ice crystal size in- 
cluded in the calculation. For example, the radiative forcing 
values at 50 /xm on the abscissa were calculated using only 
crystals in the size distribution larger than 50 /xm. It turns 
out that the solar and infrared radiative forcings are approx- 
imately equal and opposite in this cloud, resulting in a near 
zero net forcing. As the minimum size included increases, 
the magnitudes of both the solar and infrared forcings de- 
crease, with solar forcing decreasing faster, resulting in a 
slight increase in net forcing. However, the bottom line is 
that ice crystals smaller than 50 /xm contribute very little to 
the solar, infrared, or net forcing. Figure 14 shows cloud ra- 
diative heating rate profiles from the same calculations, both 
including and excluding crystals smaller than 50 /xm. As ex- 
pected the small crystals also contribute very little to radia- 
tive heating in the anvil. 

Anvil cirrus sampled on several additional days indicate 
similar general features to the cases shown above: increasing 
effective radius with decreasing height in the anvils and very 
little contribution of small crystals to SAD (i.e., extinction). 
With the exception of the active turrets sampled on 22 and 24 
July and the isolated cases discussed below, we generally do 
not see regions with elevated ice concentrations that would 
be indicative of local nucleation events within the anvils. In 
other words, the anvils sampled during TC4 seems to be con- 
trolled by dissipation and sedimentation of ice detrained from 
the deep convection, without evidence for significant regen- 
eration of ice crystals downwind from the convection. 

On at least two flights (22 July and 8 August), isolated 
occurrences of high ice concentration (>1 cm -3 ) were indi- 
cated by the 2D-S data in aged anvil cirrus far away from 
any active convection. On each of these flights there are 
a few brief periods (<lmin of flight time ~ 1.2 km dis- 
tance) with ice concentrations exceeding 1 cm -3 and some- 
times approaching 10cm _1 . Nitric oxide, carbon dioxide, 
and aerosol measurements made on the DC-8 give no indica- 
tion that these events were aircraft contrails. 

One of these high ice concentration events occurred near 
cloud top allowing us to use the ER-2 remote-sensing mea- 
surements to investigate the cloud. Figure 15 shows a time 
series of ice concentration and vertical wind speed measured 
on the DC-8 for this event. The ~ 1-3 ms 1 vertical wind 
oscillations suggest the possibility of a gravity wave in the 
region with elevated ice concentration. The corresponding 
MAS image and effective radius retrievals (for the same time 
period as shown in Fig. 15. (Fig. 16) also indicates somewhat 
of a wave structure in the cloud. However, the cloud struc- 
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Fig. 12. Ice crystal size distributions in 1-km altitude bins con- 
structed from the 8 August DC-8 and WB-57 2D-S anvil cirrus mea- 
surements. 
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Fig. 13. Solar, infrared, and net top-of-atmosphere cloud radiative 
forcings are plotted versus the minimum size included in the ice 
crystal size distributions. Incrementally excluding small crystals up 
to a larger and larger size decreases both solar and infrared forcing, 
with a slight increase in net forcing. Small ( D <50 /xm) contribute 
very little to radiative forcing. 


ture looks more like small-scale convection than purely lam- 
inar flow in a wave. The MAS effective radius retrieval indi- 
cates very small ice crystals ( r e <12 /xm) collocated with the 
high ice concentrations indicated by the 2D-S measurements. 
The 2D-S effective radii are also considerably smaller in the 
region with enhanced ice concentration than in surrounding 
regions. The likely explanation for the high ice concentra- 
tions is homogeneous freezing of aqueous aerosols result- 
ing from rapid cooling either in gravity waves or convective 
updrafts. 
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Fig. 14. Solar, infrared, and net cloud radiative heating rate profiles 
are plotted from calculations including all crystal sizes (solid) and 
excluding crystals smaller than 50 fim (dashed). The small crystals 
have negligible impact on radiative heating. 



Fig. 15. Time series of vertical wind speed measured by MMS 
(red) and 2D-S ice concentration (black) in a region with high ice 
concentrations apparently caused by convection or gravity waves. 


We emphasize here that these high ice concentration 
events are very rare in the aged anvil cirrus sampled in 
TC4. The vast majority of the TC4 2D-S ice concentration 
measurements in anvil cirrus indicated values between about 
0.01 and 0.4 cm -3 (Fig. 17). In fact values exceeding 1 cm -3 
occurred less than 1.4% of the time in anvil cirrus. As a re- 
sult, the enhanced ice concentration events probably had little 
impact on either the anvil radiative properties or persistence. 


5 Summary and Discussion 

The TC-4 microphysics dataset provides multiple lines of 
evidence suggesting that the CAS measurements of small 
crystals in anvil cirrus are strongly influenced by artifacts 
of large-crystal shattering on the instrument inlet. First, we 
compared the CAS small-particle measurements with those 
from the 2D-S instrument. The 2D-S probe is physically de- 
signed to limit shattering, and interarrival times are used to 
identify and remove shattering artifacts in the 2D-S dataset. 
The CAS small-particle concentrations are consistently 1-2 
orders of magnitude higher than the 2D-S small particle con- 
centrations. The discrepancy extends throughout the over- 
lapping size range for the two instruments (~ 10-50 /im). 

Second, we examined the relationship between the ratio 
of small-crystal concentrations indicated by CAS and 2D-S 
and the large-crystal mass. These two quantities are strongly 
correlated: as large-crystal mass increases, the CAS/2D-S 
concentration ratio increases. This dependence is consistent 
with shattering contamination in the CAS dataset. There is 
no known reason to expect that an undercounting of small- 
particles by the 2D-S probe would result in such a correlation 
with large-crystal mass. 

Fastly, we examined the correlation between small-crystal 
concentration and large-crystal mass indicated by the dif- 
ferent instruments. As anvils age, gravitational size sorting 
tends to separate large and small crystals. Hence, in aged 
anvil cirrus, the expectation is that there will be regions with 
plentiful large crystals and relatively few small crystals. In 
the DC-8 2D-S dataset such regions are apparent, and the 
correlation between small-crystal concentration and large- 
crystal mass is relatively weak. In contrast, the CAPS dataset 
indicates a strong correlation between small-crystal concen- 
tration (indicated by the CAS probe) and large-crystal mass 
(indicated by the CIP probe), as well as an absence of regions 
with abundant large crystals and relatively few small crystals. 
The correlation indicated by the CAPS data is to be expected 
if the small-crystal concentration were strongly enhanced by 
large crystal shattering. 

McFarquhar et al. (2007) reported very similar findings 
when comparing CAS measurements with the CDP measure- 
ments of anvil cirrus in the TWP-ICE campaign. Taken to- 
gether, these results provide compelling evidence indicating 
that the CAS measurements of small crystals are substan- 
tially enhanced by shattering artifacts in tropical anvil cir- 
rus. We suggest that past CAS cirrus datasets should be con- 
sidered suspect when large crystals are present. In particu- 
lar, the CRYSTAF-FACE CAS measurements in anvil cirrus 
are likely strongly enhanced artifacts, and the findings based 
on those measurements should be revisited. Shattering will 
likely have a significantly reduced effect on the CAS probe 
when large ice crystals are predominantly absent. However, 
quantifying the small-crystal artifact production from large- 
crystal shattering will be challenging given dependencies on 
a number of factors, including crystal habit, airspeed, and 
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angle of attack. As an approximate measure for quantifying 
the impact of shattering on CAS measurements as a function 
of large-crystal mass present, scatter plots and linear fits are 
presented (Fig. 6) for spurious CAS concentration, extinc- 
tion, and ice water content versus large-crystal IWC. 

We also showed a significant correlation between small- 
crystal concentration and large-crystal mass in the WB-57 
2D-S anvil cirrus measurements. We hypothesize that this 
correlation was caused by shattering artifacts from large 
crystals impacting the leading edge of the wing and subse- 
quently being swept back to the 2D-S probe mounted 2 m 
back under the wing. These artifacts cannot be identified and 
removed using interarrival times because they become spa- 
tially separated in the sheared flow under the wing. If this 
process is indeed occurring, then past microphysics measure- 
ment made with instruments at the wing-hatch location were 
presumably also affected by shattering artifacts. 

Assuming the DC-8 2D-S measurements provide reli- 
able (free from shattering artifacts) information about small 
(<50/xm) crystals, we used this dataset to evaluate the con- 
tribution of small crystals to extinction, radiative forcing, and 
radiative heating in the anvil cirrus sampled during TC-4. 
Examining a number of cloud systems ranging from active 
turrets to aged anvils, we find that small crystals contribute 
relatively little (<10%) to the total extinction. Using the 
2D-S size distributions measured throughout the depth of an 
aged anvil in radiative transfer calculations, we find that the 
small crystals have very little impact on cloud radiative forc- 
ing or radiative heating profiles in this particular cloud. 

The 2D-S small ice crystal concentrations in somewhat 
aged (~30 min to 3 h) anvil cirrus are typically between 0.01 
and 0.4cm -3 . Concentrations greater than 1 cm -1 do occur 
in the aged anvil cirrus, but they are rare (<1.4% of the aged 
anvil cirrus measurements). The measured ice concentrations 
in aged anvil cirrus are substantially lower than those mea- 
sured in convective turrets and fresh outflow (1-30 cm -3 ). 
We argue that sublimation and deposition driven by relative 
humidity variability in the anvils, along with other processes, 
likely depletes the small-crystal population rapidly. This re- 
sult may be important for evaluation of what aerosols domi- 
nate the production of ice crystals that persist in aging anvil 
cirrus. Aerosols entrained into convective updrafts in the 
middle and upper troposphere that subsequently act as sites 
for formation of droplets and ice crystals will produce hy- 
drometeors later than aerosols in the lowermost troposphere. 
As a result, the ice crystals produced by middle and upper 
tropospheric aerosols may produce the smallest ice crystals 
detrained into anvil cirrus that will be the first to sublimate 
as dry air is entrained. 

The findings reported here do not imply that small crys- 
tals don’t play a significant role in the radiative properties 
of other types of cirrus, such as anvils generated by con- 
tinental convection, midlatitude cirrus, or even other tropi- 
cal regions. Numerical simulations have suggested that ho- 
mogeneous freezing of aqueous aerosols can generate large 
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Fig. 16. Modis Airborne Simulator RGB image (left) and retrieved 
effective radius (right) for the time period on 22 July corresponding 
to Figure 15. The cloud structure suggests small-scale convection. 
The black line through the center of the image shows the DC-8 flight 
track, which was essentially directly below the ER-2. 
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Fig. 17. 2D-S ice number concentration frequency distributions 
from all TC-4 DC-8 anvil sampling periods (over 13 h of 1-Hz data). 
The convective turrets sampled on 22 and 24 July are excluded. Ice 
concentrations exceeding 1 cm -3 occur less than 1.4% of the time. 

ice concentrations in very high, cold cumulonimbus clouds 
Jensen and Ackerman (2006). Also, the uppermost tropi- 
cal tropopause is often saturated with respect to ice (Jensen 
et al., 2001), perhaps allowing small crystals to persist even 
as ambient air is entrained into the anvils. The TC-4 anvil 
cirrus measurements discussed here were relatively low and 
warm, where the ambient air is typically significantly subsat- 
urated. Knollenberg et al. (1993) used a modified 2D-C to 
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get forward-scattering spectrometer measurements of small 
crystals near the tops of very high, cold tropical cirrus. They 
reported very high concentrations (> 10 cm -3 ) of small crys- 
tals. Given the general lack of crystals larger than ~ 1 00 /i m 
in the cloud-top regions sampled, it is possible that these 
measurements were not substantially contaminated by shat- 
tering artifacts. However, it is difficult to evaluate the ac- 
curacy of these measurements since the modified 2D-C was 
used only during the one field experiment. Recent FSSP 
measurements in high, cold tropical anvil cirrus indicated ice 
concentrations in the 0. 1-1 cm -3 range (Kramer et al., 2008). 

As discussed in the introduction, accurate measurements 
of cirrus ice crystal size distributions are necessary for evalu- 
ation of cirrus radiative effects, development and evaluation 
of remote-sensing algorithms, evaluation of aerosol impacts 
on cirrus, and ultimately representation of cirrus in climate 
models. The results of this study, along with other recent 
analyses, suggest that much of the existing database of cirrus 
measurements of small crystals may be suspect because of 
the shattering problem, at least under some conditions. This 
situation motivates a need for additional airborne cirrus mea- 
surements using recently-developed probes that are designed 
to limit shattering artifacts, along with post-processing anal- 
ysis techniques to identify and remove the shattering prob- 
lem. In particular, it would be valuable to make measure- 
ments in various types of cirrus, including synoptically gen- 
erated mid-latitude cirrus, anvil cirrus produced by continen- 
tal convection, and anvil cirrus in the cold, uppermost tropi- 
cal troposphere. 
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